function d = hybridNeuron(t, y, params)
    % hybridNeuron 神经元模型的微分方程
    % 状态向量 y [x; y; z]
    % params 为参数结构体，包含字段：
    %        a, b, c, alpha, beta, delta, A, f
    %
    % 对应论文公式(5):
    
    x     = y(1);
    y_var = y(2);
    z     = y(3);
    
    % 外部激励信号
    u_s = params.A * cos(params.f * t);
    
    dxdt = (1 - params.a)*x - (1/3)*x^3 - y_var + u_s;
    dydt = params.c * (x - params.b * sin(params.beta * z) * y_var);
    dzdt = params.delta * y_var - params.alpha * z;
    
    d = [dxdt; dydt; dzdt];
    
end
